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Abstract 

We present generalizations of the spin Sutherland systems obtained earlier by Blorn 
and Langmann and by Polychronakos in two different ways: from SU(n) Yang-Mills 
theory on the cylinder and by constraining geodesic motion on the IV-fold direct product 
of SU(n ) with itself, for any N > 1. Our systems are in correspondence with the Dynkin 
diagram automorphisms of arbitrary connected and simply connected compact simple 
Lie groups. We give a finite-dimensional as well as an infinite-dimensional derivation 
and shed light on the mechanism whereby they lead to the same classical integrable 
systems. The infinite-dimensional approach, based on twisted current algebras (alias 
Yang-Mills with twisted boundary conditions), was inspired by the derivation of the 
spinless Sutherland model due to Gorsky and Nekrasov. The finite-dimensional method 
relies on Hamiltonian reduction under twisted conjugations of A-fold direct product 
groups, linking the quantum mechanics of the reduced systems to representation theory 
similarly as was explored previously in the N = 1 case. 
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1 Introduction 


The Calogero-Moser-Sutherland type many-body systems [DM and their generalizations 
continuously attract attention since they are ubiquitous in physical applications and are related 
to intriguing mathematical structures. See, for instance, the reviews MM- A powerful 
approach to these integrable systems consists in viewing them as Hamiltonian reductions 
of higher dimensional free systems. This method first came to light in the papers [8] |9]. In 
particular, Kazhdan, Kostant and Sternberg [9] obtained the classical Sutherland Hamiltonian 


U(q,p) 


Ifpl+T, 

k=l ijpj 


sin 2 (qi - q.j) 


( 1 . 1 ) 


by reducing the kinetic energy of the geodesic motion on the group SU{n) using the symme¬ 
try defined by the conjugation action of SU(n ) on itself. Their reduction required fixing the 
Noether charges of the symmetry in a very special manner. Later it turned out that more 
general choices lead to extensions of the Sutherland system by ‘spin’ degrees of freedom that 
belong to reductions of coadjoint orbits of SU(n ) by the action of the maximal toral sub¬ 
group. Spin Sutherland systems of this kind appear for every simple Lie group pnunmug. 
Moreover, not only the choice of the constraints, alias the value of the momentum map, but 
also the underlying symmetry can be generalized in several ways. For example, we obtained 
spin Sutherland systems by reducing the free motion on simple compact Lie groups utilizing 
twisted conjugations D3- These reductions were analyzed m at the quantum mechanical 
level, too, which yields a bridge from harmonic analysis to integrable systems. 

The derivation of the Calogero and Sutherland systems given in [9] represents a paradigm 
for the reduction approach to integrable systems. This framework has since been expanded 
by the introduction of several new parent systems as master integrable systems. We mention 
only the influential work of Gorsky and Nekrasov cans], where the Sutherland Hamiltonian 
was re-derived from a free Hamiltonian on an infinite-dimensional phase space built on su(n )- 
valued currents on the circle. As discussed also in M, this connects the Sutherland system 
to SU[n ) Yang-Mills theory on the cylinder. 

Our present work was motivated mainly by questions stemming from the papers of Blorn 
and Langmann [ TH1 [T9] and Polychronakos [20], where a family of generalized spin Sutherland 
systems was derived by means of two different methods. In these systems the particle posi¬ 
tions are coupled to spin variables living on N arbitrary coadjoint orbits of SU(n), and the 
Hamiltonian also involves N arbitrary scalar parameters, for any integers n > 1 and N > 1. 
Blom and Langmann obtained their systems from SU{n ) Yang-Mills theory on the cylinder, 
placing non-dynamical ‘color charges’ at N arbitrary locations on the circle. Thus their deriva¬ 
tion fits in the Gorsky-Nekrasov framework. On the other hand, Polychronakos proceeded by 
constraining the geodesic motion on the Y-fold direct product of SU{n ) with itself, using N 
arbitrary scale factors in the definition of the bi-invariant Riemannian metric. Although it was 
not stressed in [20], this amounts to Hamiltonian symmetry reduction with respect to twisted 
conjugations on the direct product group, defined with the aid of the cyclic permutation of 
the Y-factors. Based on direct comparison of the Hamiltonians, both Blom-Langmann [19] 
and Polychronakos ra pointed out that their respective systems coincide, but they did not 
provide any conceptual explanation of this remarkable fact. It is natural to ask for a better 
understanding of the mechanism behind this coincidence. 

In this paper we describe two derivations of group theoretic generalizations of the systems 
of PS, 20] and shed light on the mechanism whereby these two derivations lead to the same 
outcome. The first approach relies on symplectic reduction of geodesic motion on the Y-fold 


2 



direct product G x ■ ■ ■ x G for any compact simple Lie group G. The underlying symmetry 
is built from twisted conjugations involving the cyclic permutation of the N factors and, 
for simply laced groups, also a Dynkin diagram automorphism of G. For N — 1 these are 
the reductions analyzed in m mj. The second method utilizes a free Hamiltonian on the 
cotangent bundle of an infinite-dimensional configuration space consisting of Lie algebra valued 
(quasi-)periodic currents. The reduction is based on the natural action of the corresponding 
(twisted) loop group, where the twisting can be non-trivial in the simply laced cases. This 
second derivation fits in the framework of Gorsky-Nekrasov naunns] and Blom-Langmann 
[18, US]. In the infinite-dimensional reduction a partial gauge fixing associated with a finite¬ 
dimensional gauge slice will be exhibited that will be explicitly mapped into a gauge slice of a 
corresponding partial gauge fixing that arises in the finite-dimensional reduction. The mapping 
between these gauge slices will be seen to be a bijection that engenders an isomorphism 
of the respective reduced systems. When G = SU(n) equipped with the identity diagram 
automorphism, then the reduced systems are precisely those found in [18], 20]. 

The organization of the paper is as follows. Section 2 is devoted to group theoretic prelimi¬ 
naries regarding the relevant twisted conjugations of the direct product groups. No new results 
are claimed in this section, even though we could not find references for the statements of the 
two technical lemmas that we prove. Section 3 contains the analysis of the finite-dimensional 
Hamiltonian reduction under twisted conjugations. The structure of the reduced system is 
given by Proposition 3, which represents our first new result. The formula of the reduced 
free Hamiltonian displayed by equation (13.23P is spelled out in detail in Subsection 3.2 in the 
case of the identity diagram automorphism for any group G. Proposition 4 implies that this 
reproduces the spin Sutherland systems of [IB[ ED] for G = SU(n). In Section 4 we present 
an alternative derivation in the general case based on (twisted) current algebras. Our second 
result is that we here explain the equivalence of the reduced systems that result from the 
finite-dimensional and from the infinite-dimensional derivations. We conclude in Section 5 
by outlining the quantum mechanical analogue of the finite-dimensional classical Hamiltonian 
reduction of Section 3. There is also an appendix, where a useful technical result is described. 


2 Twisted conjugations and the monodromy matrix 

In this section we collect the group theoretic results that will be needed later. 

Let G be a connected and simply-connected compact simple Lie group. Choose an au¬ 
tomorphism 7 of G that corresponds to a diagram automorphism, denoted by 7', of the Lie 
algebra Q of G with respect to the Cartan subalgebra T < Q. If Q is not simply laced then 7 
is the identity automorphism, otherwise it can have order 1, 2 or 3; the latter occurs only for 
G = Spin(8). The group G acts on itself by the twisted conjugations G 7 that for each ij 6 G 
operate as 

C]:G^G, ( 2 - 1 ) 

The space of the corresponding G-orbits can be identified Ea with the space of the orbits 
of the so-called twisted Weyl group acting on the fixed point set T 7 of 7 in the maximal 
torus T < G associated to T. For non-trivial 7, the G-orbits in question are termed twisted 
conjugacy classes. The maximal (also called principal) orbits with respect to the G 7 action 
fill a dense open subset G' C G, and the space of these orbits can be parametrized by the 
interior of a convex polytope, T 7 C T 7 . This means that every G 7 orbit of G in G' passes 
through a unique element of the form for some /i e T 7 . Moreover, the isotropy subgroup of 
these elements equals T 7 , that is, if n £ T 7 then G 7 (e M ) = holds if and only if r) 6 T 7 . The 
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closure of T 7 is often called ‘Weyl alcove’, since it is a fundamental domain of the (twisted) 
affine Weyl group acting on T 7 . Its explicit form can be found for example in [22] (see also 

Hi)- 


Now pick a positive integer, N, and consider the direct product 

N times 


t-.= l 


Gx-xG. 


( 2 . 2 ) 


Adopting an admittedly unusual (but convenient) notation, the elements of will be encoded 
as TV-component ‘vectors’ 


~t = 9i © ■ ■ • 0 9n, 9k eG, Wk = 1,..., N. 


(2.3) 


The group operation is of course given by componentwise multiplication. We denote the Lie 
algebra of G^ as 

N times 


(2.4) 


with elements 


A — X\ © • • • © X jv, A/; G Q, Wk — 1,. .., N. (2.5) 

By combining 7 G Aut(G') with the cyclic permutation automorphism of ff we introduce 
T G Aut(Cr) by 

r (? 7i © 772 © - - - © tjn) ■= 7 (Vn) © r/1 © ■ ■ ■ © Vn-i ( 2 . 6 ) 

and let T' denote the corresponding automorphism of . We then dehne an action of ff 
itself by the T-twisted conjugations that for each if G G operate on G as follows: 


on 


C^(t) := r(if)lflf ^ = n/iVN^iVi 1 © V192V2 1 © ■ ■ • © Vn-i 9 nVn- 


(2.7) 


The T-twisted conjugations on and the 7 -twisted conjugations on G are related by the 
‘monodromy map’ M: G —> G given by 


:= gig 2 ■■■gN- 

It is easy to check the equivariance property 

M(C r lf (t)) = Cf N (M(t))- 


( 2 . 8 ) 


(2.9) 


This implies immediately that if if and if' lie on the same orbit with respect to the C r action 
of 3, then their ‘monodromy matrices’ lie on the same orbit of G under the C 7 action. The 
next lemma states that the converse is also true. 

Lemma 1. The elements if and 1 f' lie on the same orbit offf under T-twisted conjugations 
if and only if M (if) and M( g') lie on the same G-orbit under 7 - twisted conjugations. Fur¬ 
thermore, the corresponding isotropy subgroups of if inff arid M( if) in G are isomorphic. 

Proof. Suppose that if, ~ff' G G verify M(lf') = C 7 jv (M(^’)) for some g N G G. After fixing 
such an element i] N , the requirement g' = C^-(lf), which can be spelled out as 


g[ = iMgiVi , g' 2 = vid2V 2 , ■ ■ ■, 9n-i = vn-29n-iVn-i> 9n = VN-igNV N \ ( 2 - 10 ) 
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permits us to uniquely determine rj i,... ,r] N _ 1 . In fact, we simply solve the last (N — 1) of 
these equations for t)n-i, • • •, Vi, and then the first equation is already guaranteed by M(g') = 
The same calculation for shows that Miff) = Cf N (M(~f)) if and only 

if^ = C^) with 

Vn-i = 9nVn9n, Vn -2 = ( 9N-i9N)VN(gN-igN ) \ ■ ■ ■, Vi — (92 ■ ■ ■ 9n)vn(92 ■ ■ ■ 9n) \ 

( 2 - 11 ) 

which explicitly establishes the isomorphism claimed by the lemma. □ 

To continue, fix positive numbers A& (k — 1,..., N) and define an invariant scalar product 
on ~(j by 

N 

:=^A k (X k ,Y k ), (2.12) 

k =1 

where ( , ) is a multiple of the Killing form of Q, i.e., (X,Y) = — C tr(adx 0 ady) with a 
constant C > 0. For later convenience, we assume that 

y = 1. (2.13) 

Let T /T be the transpose of T' with respect to this scalar product. The defining property 

(?, r'(?)7 = (r' T (?), L) t , (2,i4) 

entails that 

T /r : Ad © • • • ® A0v-i © X N hg W © • • • © X N © ^(V)- 1 ^), (2.15) 

Ai Aat_i An 

where we used that 7 ' preserves the Killing form of Q. 

We now introduce two Abelian subalgebras fC and Q of y, which are isomorphic to T 7 : 


£~Tg HS = Vt\X 1 = X 2 = - 

• • = x N g r 7 }, 

(2.16) 

S.= {fl«ft = y. 

A k 

(Vfc = 1,... A)}. 

(2.17) 

We also define 

C := { 9> 1 9* = A 9 

Ak 

(V/c = 1,... A)}. 

(2.18) 

It is straightforward to verify that 



Q C Ker(T ,T — e _ad T) and Im(T ,T 

-e- ad ^)C/C x , Vg G T 7 , 

(2.19) 


where ‘perp’ refers to orthogonal complement with respect to the scalar product (12.12j) . The 
following strengthening of these properties will be crucial in our considerations. 

Lemma 2. For any q G T 7 , parametrizing ~f as in H2.18 1) , we have 

Ker(T ,r - e" ad T) = Q and Im(T /T - e“ ad T) = /C x . (2.20) 

Consequently, the map 

Z{q) := (r ,T -e- ad ^)| fi x: Q x ^/C x (2.21) 

is a linear bijection for any q G T 7 . 
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Proof. On account of the first relation in (I2.19p and standard linear algebraic facts, using also 
that dirn(Q) = dirn(/C), it is sufficient to prove that 


Ker(r'-e ad ^) = 1C, Vg G T 7 . 
Suppose that X G y belongs to this kernel, which means that 

r(?) = e ad Tp?). 

By taking exponential in the group Gr, we conclude that 


r(e^) = = e^e^e - ^, 


Vi G 


( 2 . 22 ) 


(2.23) 


(2.24) 


This in turn means that e ^ belongs to the C r isotropy group of . 

To finish the proof, note that M(e = e q holds for all q G T 7 . Therefore we see from (the 
proof of) Lemma 1 that if q G T 7 , then the isotropy subgroup of e^ with respect to the C T 
action is precisely the connected subgroup K of G associated with the Lie algebra /C (alias 
the diagonal embedding of T 7 into Gr). 

By combining the above, we have shown that (12.231) with q G T 7 implies that X G 1C, i.e., 
the relation ( 12 . 221 ) is valid. □ 


Remark. For q G T 7 , consider the transpose of Z(q): 

Z(q) T = (T'-e ad T)|^ : K L ->• Q ± . (2.25) 

Translated into geometric terms, the equality Im(^(g) T ) = Q 7 says that 

T e yC§(e^) = Q ± e^. (2.26) 

Here C^(e^) is the C^-orbit through e^ and on the right-hand-side we notationwise ‘pretend’ 
that we are dealing with a matrix Lie group. Therefore we have the direct sum decomposition 

T e ~t~Cj = = Qe^ ® Q ± e t = T e yQ ® 7>G7§(e^), G Q, (2.27) 

where Q is the connected subgroup of with Lie algebra Q. This is an orthogonal decompo¬ 
sition with respect to the bi-invariant Riemannian metric on G^ induced by the scalar product 
(I2.12P on . It is also worth noting that the map Z{q) is equivariant under the actions of 

K = yjiag (2.28) 

given by the corresponding restrictions of the adjoint action of G^ on y . 


3 Hamiltonian reduction based on twisted conjugations 

We are interested in reductions of the free particle moving on the group manifold G^ equipped 
with the bi-invariant Riemannian metric belonging to the scalar product (I2.12p on Q. This 
metric is clearly invariant also under the twisted conjugations (12.71) . We shall characterize 
the reductions with respect to this symmetry using the standard framework of symplectic 
geometry [23] , 
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3.1 General description for any 7 £ Aut(G) 


We begin by noting that the phase space of the geodesic motion is the cotangent bundle of 
G. We adopt the identification T*(^ ~ x y set up utilizing right-translations and taking 
as the model of by the scalar product ( , W. We also choose a coadjoint orbit of cf, 


~6 := 01 


o 


N, 


(3.1) 


where the Ok are coadjoint orbits of G. We then consider the extended phase space 

p ■= x ~6 ~ x x ~6 = {(t,-},-£)}■ (3.2) 

Like in Section 2, we have = J\ © ■ • • © Jn and £ = £1 © ■ ■ • © £at. The phase space P 

carries the symplectic form 

N 

12 = d(3, dlflf~ l )-£ + u with u> = ^^ujk, (3.3) 

k =1 


where Uk is the natural symplectic form of Ok- The orbital part u ; implies the Poisson brackets 

«T, (t, 7>t) = (f, 7,7> ? , ( 3 . 4 ) 

for any constants X, Y E~y . The other iron-trivial Poisson brackets between the fundamental 
variables can be displayed as 

{y,<7,7) t } = 7y, {<7, W, A Ay} = 7,7,7b?. ( 3 . 5 ) 

Here, we have Xlf = X\g\ © • • • © X^g^. The Hamiltonian H of the free particle, trivially 
extended to P, reads 

= (3.6) 

The free flow generated by H through the initial value (if, J , £ ) is 

©w, d(«)7 w) = (©y. 7, f )■ ( 3 . 7 ) 

The twisted conjugations lift to the Hamiltonian action C r of C on P: 

C ^(7©?) = (C%(y), Adri^f"?), Ad^(7))- (3.8) 


This action is generated by the equivariant momentum map H/: P —> Q given by 

t) = T ,T 0) - AA-f-10) + ?. (3-9) 

Indeed, the Hamiltonian vector field of the function Ty = (T,^)y coincides with the in¬ 
finitesimal generator of the G-action (13.81) . 

According to the standard shifting trick [23], we can represent any reduction of the free 
motion with respect to the twisted conjugation symmetry as symplectic reduction of the 
Hamiltonian system (P, 12, H) at the zero value of the momentum map; in association with 
some . The resulting reduced Hamiltonian system lives on the space of t^-orbits 

Pred = Pp=o/ (3.10) 
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which inherits a reduced symplectic structure, h2 re( j, and Hamiltonian, H ie d, from and H on 
P. To be more precise, one should note that P ie d is not always a smooth manifold, but is in 
general a so-called stratified symplectic space, i.e., a disjoint union of symplectic manifolds of 
various dimensions [23] . The detailed structure depends on the choice of the orbit S. We 
here will not dwell on this issue, but shall instead focus on a dense open subset P' of the phase 
space P whose reduction can be characterized rather directly. 

Remember that every element of (it can be transformed by twisted conjugation into a 
unique element of the form e~^, parametrized according to (I2.18j) by q from the closure of the 
generalized Weyl alcove T 7 . From now on we restrict our attention to the subset P' consisting 
of triples (if, J , £ ), where the twisted conjugates of if can be parametrized by elements q 
of the open alcove T 7 . It is then clear that every f^-orbit in P ^ =0 intersects the subset S C P 
defined by 

s ■= {(e* | * (e* = get 7 }. (3.11) 

Observe also from Lemma 1 that any G-'-orbit in P^=o intersects S in an orbit of the subgroup 
K < G displayed in (I2.28I) . In other words, S is the ‘gauge slice’ of a partial gauge fixing of 
the (^-action on P^ = 0 , for which the ‘residual gauge transformations’ belong precisely to the 
subgroup K. In this way we obtain the identification 

PL d ~ PL o/l - S/K. (3.12) 

Consequently, the G -invariant smooth functions on P^ = 0 , which descend to smooth functions 
on P' ed , are equivalent to the /T-invariant smooth functions on S. The Poisson brackets of 
these functions can be determined with the aid of the restriction of the symplectic form H to 
S. 

In order to make the above concrete, we proceed to solve the restriction of the momentum 
map constraint on S. This amounts to the following requirement: 

(T /T - e -ad T) (j) + ~i = 0. (3.13) 

To solve this equation we use the decompositions 

/C + ZC^cf = Q + Q ± (3.14) 


and write accordingly 


= "?/c + T/c-l and 


—J — ^ Q + ^ 


Q 1 


(3.15) 


Here, recall the definitions of /C and Q from equations (12.16(1 and (12.171) . Referring to Lemma 
2, the /C-component of the momentum map constraint (I3.13p enforces that £ £ = 0, while its 
/W-component is completely solved by 




/C- 1 -; 


(3.16) 


using the operator Z(q) dehned in (12.211) . The element q £ T 7 can be chosen arbitrarily, and 
J q is also a free variable, which we parametrize as 


J Q=P = T~ 


P_ 

A2 


p_ 

Am’ 


p g r 7 . 


(3.17) 














Since the variables q,p and £ uniquely label the points of S, we obtain the identification 
of manifolds 

S~T 7 xT 1 x0n /C- 1 ) = {{q,p, £ K x)}. (3.18) 

According to the general principles of Hamiltonian reduction, the reduced system on P r ' ed is 
determined by the pull-backs of the symplectic form fl and the Hamiltonian H to S. Let us 
denote these pull-backs by fig and Hg. In terms of the model (13.1811 we obtain 

fig = d(p, dq) + (3.19) 

Here ^-^ nlc ± is the contribution of the orbital form c a from (13.31) . The first term is the natural 
symplectic form of the cotangent bundle 

T*f 7 ~ fl X T 7 = {(g,p)}. (3.20) 

Since the residual gauge transformations by K = T diag act only on £ the reduced phase 
space can be decomposed according to 

P r ' ed ~ T*f y x (3 n K l )/K. (3.21) 


This is a symplectic identification in the following sense. The natural action of the group 
K < G on is generated by the momentum map £ n- £ A: , and 

5 red = (3n /C ± )//l (3.22) 


is the reduced orbit obtained by setting this momentum map to zero. We denote the (in general 
stratified) symplectic structure inherited by C^ re d as u> ie d. This is carried by the second factor 
in the Cartesian product (13.211) . Regarding the reduced Hamiltonian of the free motion, we 
already mentioned that it is determined by Hg. By means of the model (13.181) , Hg takes the 
form 

Hg = ~(p,p) + ~{Z{q) 1 (£/c- l ),Z(q) 1 (£/c-l))a- (3.23) 

The first term is the kinetic energy of the free motion on T 7 with respect to the flat metric, 
while the second (also K -invariant) term describes the interaction of the spatial degrees of 
freedom q,p with the orbital variables given by the K -orbits in (13.221) . The outcome of the 
foregoing discussion can be summarized as follows. 


Proposition 3. A model of the open subset of the reduced phase space / 1 3.1(A) associated with 
generic monodromy matrices parametrized by q e T 7 , is provided by P*T 7 x d red • 

Here T *T 7 carries its Darboux form and ~3 re d jK 2 A) is equipped with its own (stratified) 
symplectic structure arising by reduction via the K-action. The reduced free Hamiltonian on 
this subset P r ' ed = S/K is encoded by the K-invariant function Hg \3.2S 1) . 

We finish this subsection with a few comments. First, we remark that the Poisson bracket, 
{ , } re d, that is induced on the smooth functions on P r ' ed is encoded by 


{(q,X), (p,Y)} red = (X,Y), VX,y e r\ (3.24) 

together with the restriction of the Poisson brackets of the /T-invariant functions on to 
D IC L , which can be determined by means of the Lie-Poisson brackets (13.41) . 

Second, we recall that the N — 1 special case of our classical Hamiltonian reduction 
was previously studied in [13] and the corresponding quantum Hamiltonian reduction was 
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investigated in [TJj. As will be discussed in Section 5, the quantum Hamiltonian reduction 
can be described similarly in the general case. To facilitate the comparison with these earlier 
works, we now recast the inverse of Z(q) ( 12 . 211 ) in the form 

Z(q)~ l = — e ad T o ( U(q) T )-\ (3.25) 

where the linear operator U(q) : /C x —* reads 

U{q)= (id^-e- ad ^or')| K x. (3.26) 

Since e ad T preserves the scalar product on y , and inversion and taking transpose commute, 
we can rewrite the Hamiltonian H$ (I3.23P as 

Hs — ^(p,p) + k^-i 1 ('Ck:- l ))t- (3.27) 

For N = 1, when Q = /C, it was shown in [13] that the extension of U(q)^ 1 by zero on Q is 
a solution of the modified classical dynamical Yang-Baxter equation. There it was important 
that the automorphism y' preserves the scalar product on the underlying Lie algebra. This 
holds for our T' if and only if Ai = A 2 = . . . = Atv- It would be interesting to know if U(q) is 
still related to a classical dynamical r-matrix for any N with arbitrary parameters A . 

Third, it is proper to mention here that Hochgerner [[12] studied reductions of cotangent 
bundles under the assumption that a single isotropy type occurs for the underlying group 
action on the configuration space. Our Proposition 3 can be obtained as a special case of his 
results. Since this would necessitate going into several technicalities, we have preferred to give 
a simple direct derivation. 

Finally, note that the reduced system governed by the Hamiltonian (I3.27j ) is exactly solvable 
since its solutions can be obtained by the projection method applied to the free flow (I3.7[) , 
which involves only algebraic operations. There exist general arguments Gam on! that also 
indicate Liouville integrability on the full reduced phase space P re d- Indeed, one can easily 
construct many conserved quantities in involution as follows. For each fixed real parameter 
u 7 ^ 0 , define the -valued function 0 U on the unreduced phase space by 

+ (3.28) 

u 

Taking a basis {V a } of y , the Poisson brackets of the components 0“ := (0 U , V a )-g satisfy 

{€, £} = St (—€ + —k) (3.29) 

\U — V V — U J 


for any u ^ v, where [H a , V b ] = f^ b V c and summation over the index c is understood. It 
immediately follows from (13.29)1 that for any two (^-invariant smooth functions hi and h 2 on 
y and for any parameters u and v we have 


{hi o 0 U , h 2 o cj) v } = 0. (3.30) 

This descends to a Poisson commuting family on the reduced phase space. As generators, it is 
useful to take invariant homogeneous polynomials A and extract the coefficients of the powers 
of u^ 1 from hi o (j> u . Because (13.30(1 holds for any u and v, the so obtained coefficient functions 
are also in involution. The w-independent term coming from h(X) = gives the 

main Hamiltonian after reduction, and it is natural to expect that the family just delineated 
is generically sufficient to ensure its Liouville integrability. 
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3.2 The examples associated with 7 = id^ 

Our goal now is to provide an explicit formula for the reduced Hamiltonian Hg (I3.27|) in the 
case of the identity automorphism 7 of an arbitrary connected and simply-connected compact 
simple Lie group G. During the calculations it proves to be handy to realize the real Lie 
algebra Q as a compact real form of its complcxihcation Q c . Since the complexification T c of 
T is a Cartan subalgebra of Q c , the pair ( Q c , 7” c ) uniquely determines a reduced root system 
<f>. Choose a polarization $ = <f> + U <f>_ and let n C <f> + denote the corresponding set of simple 
roots, whose number is the rank of G, r. Also, for each root 99 we select an appropriate root 
vector X lf . On the root vectors we may and shall impose the following two conditions: 

(a) The root vectors are normalized by the conditions 

k(X v ,X_ v ) = 1 (v?e$+), (3.31) 

where k denotes a positive multiple of the Killing form of Q c . 

(b) The real Lie algebra Q decomposes as 

Q = T© (® ve * + R Y v ) © (®„ e * + RZ v ) , (3.32) 

where for each 99 £ <L + we define 

Y v = i Xv and Z f = Xv ~J^- v . (3.33) 


Take an arbitrary set of vectors {2}}£ =1 C T satisfying n(Tj, T^) = —Sj.k (1 < j,k < r). 
Clearly the family of vectors Tj (1 < j < r ) and Y :fl Z^ (99 £ 4> + ) gives an orthonormal basis 
in Q with respect to the Ad-invariant Euclidean scalar product obtained by restricting the 
complex bilinear form ( , ) — —k onto Q. Obviously, this convention agrees with the one used 
in (12.12[) and the definition of ( , )-& extends by linearity to the complexification of £/". 

Recalling the form of the Hamiltonian Hg (13.271) . first we wish to examine the action of 
the linear operator U(q) T V((q) (q £ T, an open Weyl alcove) on a convenient basis of /CA For 
this reason, for all V £ Q c and 1 < I < N we define the vector V^ 7 £ with components 
(t’)j = S'jV (1 < J < N). To clarify our notations, we remark that for A^ (12.51) we have 
(X)i = Xj. The family of vectors 

'f/ {1 < j < r, 1 < I < N) 

gives a basis in the complex Lie algebra y c . 
the complexification of the linear operator 


and 


V 


(99 £ 4>, 1 < / < N) (3.34) 

Now we record that for each q £ T the action of 


1 (g) := (id^ - e“ ad A o r) T o (id^ - e~^ o T') : 


(3.35) 


on the basis (13.341) reads 


N 


N 


(3.36) 


l(q)f/ = T-M lj ( 0)f/ and Z{q)^ J = V -M 7 >J (^(g))^ 7 , 

Xl i~i 

where M(ix) is the N x N complex matrix with entries 
Mi t j(ix) = (A/ + Aj+i) 61 j — Xie~ 1 Xl/Xl 5i-i y j — \r+ie lx ^ Xl+1 <A+i.j (1 < /, J < N), (3.37) 
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depending on the real parameter igR. Here and below, the capital indices are understood 
modulo N. The factor i occurs since ip(q) G iR for each q G T and <p G <3>. It is also 
worth mentioning that the function M(ix) G C NxN is 27T-periodic with symmetry 

properties 

Mjjiix) = M ItJ (-ix) = M ItJ {ix). (3.38) 

On account of the formula (13.271) and the relationship 

{UbfUb))- 1 = {T{q)\ K x)~\ (3.39) 


we need control over the inverse of M(ix). One can verify that M(ix) is invertible if and only 
if x G R \ 27 tZ. Moreover, upon setting 




1 

Afc 


and bi t j = bj — bj 


(1 <I,J<N), 


(3.40) 


for all x G R \ 27 tZ the entries of the inverse matrix V(ix) = M(ix) 1 


can be checked to be 


Vi A H) = e~ ibl ’ jX 


1 1 bj j , . . 

4 5^a72) + 1 T cot(:c/2) 



(1 <I,J<N). (3.41) 


This matrix inversion was treated in detail in [20]. We note that p(q) G i(R \ 27 tZ) holds for 
every q from the open Weyl alcove T and tp G < h. 

Turning to the study of the degenerate matrix M(0), let us endow the linear space of the 
complex column vectors <C N with the inner product 


C^xC^ (*, y) ^ x*A y G C 


(3.42) 


induced by the positive definite diagonal matrix 

A = diag(Ai,...,Ajv). 


(3.43) 


Thinking of the N x N complex matrices as linear operators acting from the left on the inner 
product space C N , from the formulae (13.371) and (13.411) we can easily deduce the following: 

(i) The operator A _1 M(0) is self-adjoint with respect to the inner product (13.421) . the kernel 
J\f = ker(A -1 M(0)) is spanned by the column vector n = [1 ,..., 1]*, and the orthogonal 
complement A l' ± of A f is an invariant subspace. 

(ii) The restriction (A _1 M(0))|^/'x: A/" 1 —> AT 1 is invertible and for its inverse we have 

((A- 1 M(0))|_^x) _1 = (vr^Avr^x, (3.44) 

where the N x N matrix 7r± denotes the orthogonal projection onto the subspace A f ± , 
whereas V is the N x N matrix with entries 

V^j = {brj)2 ~ |6/ ’ j| (1 <VJ<N). (3.45) 


The operator defined by (13.441) encodes the action of the inverse (13.39)1 on AW D T■ It is 
reassuring to see that the outcome of the elementary algebraic manipulations is consistent 
with the statements in (12.19)) and Lemma 2. 
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Finally, let us observe that due to the relations (I3.39P and (I3.36|) we now fully control the 
linear operator (U(q) T U (g)) _1 ■ This implies the explicit formula of H$ given below. 

Proposition 4. For all N > 1, at each point (q,p, £ £±) £ T x T x ~3 re d, the reduced 
Hamiltonian Hs \3.2ty associated with the trivial automorphism 7 = idc takes the form 


r N 

Hs = 2 (PiP) + 9 ■> £ £ /c x )- 


j =1 i,j =1 


AT 


(3.46) 


Hpj(ip{q)){X lf I , £ K ,- L )^(X~ lf J , £ jc 1 - ) 




I / SS # I,J=1 


with V and V defined in \3.fl\) and \3.fdd , respectively. 


The above Hamiltonian Hs can be interpreted in terms of ‘particles’ with positions defined 
by the components of q that interact with each other as well as with the ‘spin variables’ rep¬ 
resented by £ . It is a system of ‘Sutherland type’ since the interaction exhibits trigonometric 
dependence on the coordinates, which means that the ‘particles’ move on the circle. For 
Q = su(r + 1), the explicit formula (13.461) reproduces precisely the integrable spin Sutherland 
Hamiltonian obtained previously by Blom and Langmann [TH1 [19] and by Polychronakos [2D] 
by means of different methods, and generalizes their system for arbitrary simple Lie algebras. 


4 Alternative derivation from twisted current algebra 

It is remarkable that certain many-body systems can be obtained by alternative reduction 
procedures starting from finite-dimensional as well as from infinite-dimensional ‘free systems’. 
The main ideas in this regard go back to Gorsky and Nekrasov mm, who first interpreted 
the standard trigonometric Sutherland system from an infinite-dimensional standpoint. Their 
derivation had been exposed in [5] [25], too. In effect, it also formed the basis of the derivation 
of the Q = su{r + 1 ) special case of H$ (13.461) given by Blom and Langmann [181 US]- We 
thought it worthwhile to develop the Gorsky-Nekrasov method in our general case. Below we 
re-derive the same systems as those obtained in Section 3.1, but we will not be fully rigorous 
as we shall neglect topological details concerning the infinite-dimensional spaces involved. 

We start with the infinite-dimensional phase space P consisting of triples of (/-valued func¬ 
tions (A(x), E(x),(( x )) 011 the rea l li ne , subject to the following quasi-periodicity condition 
for C = A, E, 0 

C{x + 1) = r'(C(x)) with r := y” 1 . (4.1) 

We assume that £(x) is actually a generalized function of the form 

- x k ), (4.2) 

fcez 

where 

0 < xi < X 2 < ■ ■ ■ < xn < 1 (4.3) 

and 

Xk+N = Xk + 1, Ck+N = t'(C k), Vk G Z. (4.4) 

The normalization is such that J 0 ' S(x — Xk)dx = 1 for k = 1,..., N. Moreover, we assume 
that (k belongs to an adjoint orbit 0' k of G in Q , identified with Q* via the scalar product ( , ). 
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Regarding the quasi-periodic fy-valued functions A and E, we suppose that they are smooth 
on each open interval (xk,Xk+i), and they as well as their derivatives of all orders have finite 
one-sided limits at every point Xk- For reasons that will become clear shortly, we permit both 
A, E and all their derivatives to have jumps at the distinguished points Xk- 

By using smooth ^-valued quasi-periodic test functions V(x) and W{x), the non-vanishing 
Poisson brackets on P are specified as 


{A(x),V(x))dx, (E{x),W{x))dx\ = / {V(x),W(x))dx 


and 

{((k,V k ),(<k,W k )} = (( k ,[V k ,W k }), 

where 14 : = V(xk), B4 := W{xk). The corresponding symplectic form fl on P reads 

N 


Q = dd + hj' with u' = 


(4.5) 

(4.6) 

(4.7) 


k=1 


where oj' k is the Kirillov-Kostant-Souriau form on the orbit 0' k and 

(E(x),dA(x))dx. (4.8) 

We equip the phase space P with the free Hamiltonian 

(E(x), E(x))dx (4.9) 

and a Hamiltonian action of a suitable twisted loop group, G. Namely, an element g G G, 
which is a G-valued function on the line subject to g{x + 1) = r(g(x)), acts on the triple 
(H, E, () according to 

A(x) i y g(x)A(x)g~'(x) - g'fxjg^fx), E(x) gfijE^g -1 ^), ( t i-> g(x t )( lt g'\x i ). 

(4,10) 

We suppose that g and all its derivatives are smooth on each open interval (xk, Xk+i) and have 
finite one-sided limits similarly to A and E. An important difference is that g is required to 
be globally continuous on R (including at the points Xk)- These assumptions guarantee that 
the formulae in (14. lOj) yield a well-defined action of G on P. 

The action (14. 10|) admits the equivariant momentum map 

§(x) := [A(x), E(x)] + E\x) + ((x) (4.11) 

that takes its values in the distributional dual of the space of fy-valued smooth quasi-periodic 
functions. In particular, the derivative of E(x) is understood in the distribution sense, i.e., 
the value TfR] of the functional T on the test function V is 

rl N 

T[H]= / ((V(x),[A(x),E(x)}) - (V'(x),E(x)))dx + J2(V(xk),(k). (4.12) 

Jo k=i 
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By using this, one can verify that 'I' generates the infinitesimal action according to 

{A(x)MV]} = [V(x),A(x)]-V'(x), 

{E(x),y[V]} = [V(x),E(x)}, (4.13) 

{( k ^[V}} = [V(x k )Xk\. 

The reduction of our concern is defined by imposing the momentum map constraint 

T = 0 (4.14) 

and then factorizing by G , as usual. An important fact used in the subsequent analysis is that 
for every A that appears in P there exists an element g e G that transforms it into a constant 
function on M. The constant can be taken from the closure of the open ‘twisted Weyl alcove’ 
T 7 that parametrizes the twisted conjugacy classes in G given by the orbits of the C 1 -action 
m- For convenience, we present a proof of this fact in Appendix A. 

Similarly as we worked in Section 3, we restrict attention to the subset P' consisting of 
such triples for which A(x) can be gauge transformed to a constant from the interior of the 
alcove, denoted as 

X e f 7 . (4.15) 

Mimicking Section 3.1, we consider the ‘gauge slice’ S' of a partial gauge fixing in P'\ 

S := { (A, E, C) | T = 0, A(x) = X 6 T 7 }. (4.16) 

Clearly the momentum map constraint T = 0 requires that the function E(x) takes the form 

E(x) = e ~G~ x A^x E + f or Xk<x< Xfc+1 , (4.17) 

where the only objects so far undetermined are the constants 

E* e g. (4.18) 

On account of the quasi-periodicity, we have E k+N = r'{E k ) and we also define 

E~ = e -Gk- x k -i)^ xE +_^ yj, e (4. 19 ) 

Since the function E(x) has the jump (E k — E k ) at Xk, we hnd that when applied to the 
elements of S the constraint T = 0 translates into the following equations: 

( k + (E+-E k ) = 0, k = 1..., N. (4.20) 

Now, our key observation is that the correspondence 

Q = X, t- = x k ~ %k- i, A k ^k = Ck, A kh = Et- 1, k = 1,..., N, (4.21) 

*k 

allows to reformulate (I4.20P as the following system of equations: 

6 Jk+i - e~ ad ^ (J fc ) = 0, k = 1 ..., N, (4.22) 

where Jat + i = r'(Ji) and Aat + i = Ai. The crux is that this system of equations is nothing 
but the componentwise form of the relation (13. 13 j) . which represents the momentum map 
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constraint 'I' = 0 applied on the gauge slice S in the finite-dimensional symplectic reduction 
of Section 3.1. Since the residual gauge transformations are in both cases given by the action 
of the group K (12.281) . the relation (14.211) induces a one-to-one correspondence between the 
reduced phase spaces coming from the infinite-dimensional and finite-dimensional reductions: 

Ad = PkJG ^ S/K ~ S/K ~ fi, 0 /G ~ p; ed . (4.23) 

It can be readily checked that the reduced symplectic structures and also the reduced free 
Hamiltonians are converted into each other upon this correspondence. The latter statement 
follows immediately since on S we have 

„i N N 

/ (E(x), E(x))dx = y^(x fc - x k -i)(E£_ 1 , E^_f) = ("/, (4.24) 

k= l k = i 

As for the symplectic form, let us first notice that the mapping A*.£*. ^ Cfc converts the 
(unreduced) Poisson brackets in (13.41) into those in 04.61) . In other words, (Ok,u>k) and (0' k ,uj' k ) 
represent the same coadjoint orbit of G in Q*, identihed with Q via the different scalar products 
A k{ , ) and ( , ). Consequently, it is enough to focus on the 1-form 6 (14.81) . Regarding this, 
we have 

f l N 

°\s= / (E(x),dq)dx = ^(x k - x k -i)(E+_ v dq) = (f = (7 = (p,dq). 

Jo k =i 

(4 - 25) 

The second equality holds since only the T-part of E(x) contributes, which is constant on 
each interval (xk~i,Xk). We have taken into account the correspondence (14.211) together with 
equations (13.151) . (13.171) and (12.131) . We conclude that the restriction of 11 on S is mapped into 
the restriction of Q on S, which implies that the respective reduced symplectic structures are 
converted into each other, as claimed. 

The infinite-dimensional phase space that we started with emerges from Yang-Mills the¬ 
ory in (l+l)-dimensions with quasi-periodic boundary condition and external non-dynamical 
charges located at the points ,xy. Indeed, A = A\ is the spatial component of the Yang-Mills 
potential A M (/i = 0,1) in the gauge where A 0 has been set to zero, and E is its canonical 
conjugate (the ‘color-electric’ held). Because of this interesting physical context [ 151 IT8l 151 ITT] 
it would be desirable to place the above sketched derivation on completely rigorous mathe¬ 
matical ground. We plan to deal with this issue, together with the isomorphism of the full 
reduced phase spaces Pq,=o/G and P^ =0 /G, elsewhere. 

To finish, note from (14.101) that A(x) transforms in the same way as (at fixed level) the 
elements of the dual of a centrally extended twisted loop algebra based on Q. In this section 
we followed the current algebraic derivation of the Sutherland system as presented in the book 
[25], but (unlike there) we did not assume that A{x) varies in the smooth dual. Although such 
framework may appear advantageous for the symplectic reduction itself, it is problematic since 
the flow of the ‘free Hamiltonian’ H (14.91) starting at S at t — 0 leaves the space of smooth 
fields A, because E must have jumps at the points Xk- In fact, our assumptions on A, E and 
g were designed to ensure that the flow of the Hamiltonian //, 

(A(x, t), E(x, t), Ci (t), ..., GvOO) = {A{x, 0 ) + tE(x, 0 ), E(x, 0 ), (i( 0 ),..., Cjv(O)), (4.26) 
stays in the unreduced phase space P on which H is well-defined and enjoys G'-symmetry. 
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5 Conclusion 


In this paper we presented two derivations of novel spin Sutherland systems that in special 
cases were studied earlier damn Eo], and described the mechanism whereby the two deriva¬ 
tions always yield the same result. The main virtues of the infinite-dimensional derivation are 
its connection to Yang-Mills theory and that it can be extended for obtaining elliptic general¬ 
izations, similarly to the derivation of the standard elliptic Calogero system from the current 
algebra on the torus ESZl SI Id- The finite-dimensional derivation has different advantages. For 
example, it permits the construction of classical solutions by a purely algebraic projection algo¬ 
rithm. Its main advantage is that, as we outline next, the corresponding quantum-Hamiltonian 
reduction is also in the range of well-understood group-theoretic methods. 

The quantum mechanical analogue of the unreduced classical system (P, Q, H) is defined 
by the Hamilton operator 

H = A 3®idvv, (5- 1 ) 


where A-g is the Laplace-Beltrami operator associated with the Riemannian metric on G 
that corresponds to the invariant scalar product ( , W. Now this operator acts on (its usual 

domain in) the Hilbert space L 2 0 , d^)®V^. The measure dfjL-g comes from the metric and 

the vector space Ry carries a highest weight representation p-p of the direct product group 7*. 
This representation is the ‘exterior tensor product’ of representations of the iV-factors, which 
are quantum counterparts of the constituent orbits Ok of O (13. ip . i.e., 


Rp — V Vl Kl • • • KI V UN . (5-2) 

Any group element if G is represented by a unitary operator U^f on the Hilbert space, 
operating on a Rp-valued function F according to 

F^U t F = p^{lf) ofoC^. (5.3) 

The reduced Hilbert space is provided by the subspace of G^-singlets, 

(L 2 0,d^)®V-jf. (5.4) 

This is mapped to itself by Ft, and the corresponding restriction defines the reduced quantum 
Hamiltonian. By arguments similar to those in na, it is not difficult to show that the reduced 
quantum Hamiltonian is unitarily equivalent to the operator H IC d given in (15.611 below that 
acts on (a suitable dense domain inside) the Hilbert space 

L 2 {T\diir,)®V$, (5.5) 

where R^ C Rp is formed by the vectors fixed by the subgroup K (12.2811 and the measure on 

T 7 is defined by the Euclidean scalar product ( , ). To describe the operator H re d , consider 
dual bases T a and T b of ((T a , T b )-g = 5 b ) and introduce l a ’ b (q ) = ( T a ,l(q)T b )-g using l(q) 
in (I3.35p . Then one obtains 

1 cP 1 1 —> 

= -5 E f + 2 Y (5.6) 
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where the g; are coordinates on T 7 with respect to an orthonormal basis and Cq^ is a constant. 
Observe that apart from this constant (which arises from a similarity transformation [H]), the 
outcome of the quantum Hamiltonian reduction is obtained from the classical Hamiltonian 
(13.23[) via the naive quantization of the kinetic energy and the replacement of (T a , £ )-^> by 
p-p{T a ) together with restriction to the K -invariant subspace of V-j annihilated by the quan¬ 
tum analogues of the classical constraints enforcing £ £ = 0. The constant was calculate^ in 
m for N — 1 and any 7, and now we also calculated it for the systems detailed in Section 
3.2, i.e., for 7 = id^ and arbitrary N. It turned out that the constant is actually independent 

of N and of A and it reads Cq’£ = -~^k(8, 8) with the ‘Weyl vector’ 8 — \ a - 

Finally, let us sketch how representation theory can be used to determine, in principle, the 
spectrum of the reduced Hamiltonian if re d- The clue is that the spectrum of Ayt is known from 

the Peter-Weyl theorem. The theorem says that the representation of x on L 2 0 , dp,p) 
that comes from left- and right-multiplications decomposes as 


l 0,dp-g) — (Vp IE Fy) , 


(5.7) 


where all ‘components’ run over the dominant integral weights of G and A* is composed from 
the highest weights A£, denoting the highest weight of the contragredient of the representation 
V\ k . The Laplace-Beltrami operator is constant on each subspace Vp E F^, taking the value 


N 

C?(t) = -'%2 + 25, A fc ). (5.8) 

k= 1 

It is easy to see from this that the space of singlets in L 2 (G , dpp) © Vp with respect to 

the representation (15.31) of G can be presented as an infinite orthogonal direct sum of the 
finite-dimensional spaces 


(Fv* © F Al © V Ul f K • • • K (F a ^ © V Kn _ x © V Vn _ 1 ) g K (Fv* 0 y © F Ajv © V Un ) g . (5.9) 

The sum is over those vectors A for which the spaces of the G-singlets in all the N 3-fold tensor 
products in ( 15 . 91 ) are non-trivial. The composition of A^ with 7' makes sense since 7' acts on 
the Cartan subalgebra on which the weights are defined. The non-trivial finite dimensional 
spaces of the form ( 15 . 91 ) are eigensubspaces of the original Hamiltonian //, with eigenvalue 
— AC 2 a (A). It is an archetypical Clebsch-Gordan problem to determine these spaces. This 
infinite collection of Unite-dimensional linear algebraic problems is equivalent to the problem 
of diagonalization for the generalized spin Sutherland Hamiltonian H ie d ( 15 .6 |) . 
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1 A notational difference is that in mm the normalized Killing form k was denoted by ( , }. 
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A Appendix on gauge transformations 

The argument presented below is adapted from the description of the coadjoint orbits of the 
twisted affine Lie algebras [26] E], which generalizes earlier results on the untwisted case. 

Take a fy-valued function A on the real line R that satisfies the assumptions detailed in 
Section 4. In particular, it is quasi-periodic in the sense that A(x + 1) = r'(A(x)), where r' 
and r are automorphisms of Q and G induced from an automorphism of the Dynkin diagram. 
(The trivial automorphism is of course a special case.) Then consider the following differential 
equation for a G -valued function y A , 


Va(x) = Va(x)A(x), (A.l) 

where at the distinguished points Xk (see Section 4) this is understood in the sense of one-sided 
limits. We further impose the initial condition 

Va( 0) = e G G (A.2) 

and also require continuity of yA at all points x\~. It is easily seen that there exists a unique 
‘fundamental solution’ yA that meets these requirements. This boils down to well-known exis¬ 
tence and uniqueness statements on each closed interval [x fc ,x fe+1 ], and yA can be constructed 
by gluing the solutions on each closed interval by continuity. The solution yA can be writ¬ 
ten explicitly as an ordered product-integral (‘Wilson line’). It follows that yA verifies the 
conditions used in our definition of the group G in Section 4, except for the quasi-periodicity 
condition. Indeed, the quasi-periodicity of A and the uniqueness of y A imply the relation 

y A {x + 1) = y A (l)r(y A (x)), \/x G R. (A.3) 

Recall that there exists an element w G G and a unique y £ cl(7” 7 ) for which 

e x = CZ(y A ( !))• (A-4) 

Here, cl(T 7 ) is the closure of T 7 introduced at the beginning of Section 2. The automorphisms 
7 and r are inverses of each other (actually equal except possibly for G = Spin(8, R)). Picking 
w and Xi we now define the G-valued function gA on R by 

g A {x) = e _xx 7 (w)y A (x). (A.5) 

It is trivial to check that gA satisfies the identity 

g A {x)A(x)g2 1 (x) - g' A (x)g2 1 (x) = y, Vx G R, (A.6) 

and is correctly quasi-periodic 

9a(x + 1) = r(g A {x)). (A. 7) 

In conclusion, we have shown that any function A subject to our conditions defining the 
phase space P can be gauge transformed into a constant y from an arbitrarily chosen ‘twisted 
Weyl alcove’ by a gauge transformation defined by an element of the loop group G. 
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